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I. INTRODUCTION 

Frustrated spin systems have been continuously ex- 
plored in the last years. Frustration is considered a key 
ingredient to induce unconventional magnetic orders or 
even disorder, including spin-liquid states and exotic ex- 
citations. In one-dimensional and quasi-one-dimensional 
models, quantum antiferromagnets show many fascinat- 
ing magnetic properties at low temperatures which con- 
tinue to attract an intense theoretical activity. As 
representative of geometrically frustrated homogeneous 
spin chains, one can consider the antiferromagnetic spin 
S = 1/2 zig-zag chain (for which compounds such as 
CuGeOzr, LiViO^ or SrCuO^ are almost ideal proto- 
types) and three-leg antiferromagnetic spin tubes (real- 
ized in \{CuCl2tachH)^Gl]Cl^r-) . The chemistry of these 
compounds enables the synthesis of single crystals much 
larger than the previously observed organic analogs and, 
consequently, the achievement of new and more precise 
experimental studies. 

In this context, both experimental and theoretical in- 
terest on magnetoelastic chains was triggered by the 
discovery of the spin-Peierls transition in CuGeO^ at 
zero magnetization. This transition is an instability due 
to magnetoelastic effects which is characterized (below 
a critical temperature T$p) by the opening of a spin 
gap and the appearance of a dimerizcd lattice distor- 
tion at M = 0, with the consequent modulation in 
spin exchanges. Thus two related issues play together: 
the lattice distortion represents a cost in elastic energy, 
while the spin exchange modulation modifies the mag- 
netic spectrum. 

A similar phenomenon can be analyzed in magnetoe- 
lastic systems at non zero magnetization, appearing as 
most interesting the systems exhibiting magnetization 
plateaux. Moreover, it has been shown in Ref. @ that a 
spin-phonon interaction in zig-zag chains explains a spin 
gap opening as well as the presence of non-zero magneti- 
zation plateaux at low frustration, where they are indeed 
absent in the case of non-elastic chains. Such plateaux 
are due to a mechanism of commcnsurability between 
lattice distortions and spin modulation. 



Regarding non-elastic zig-zag chains at M = 1/3 mag- 
netization plateaux, it was recently shown 7 ^ that small 
modulations of exchange couplings with period three on 
top of a homogeneous zig-zag chain can drive a magnetic 
transition from a three-fold degenerate ground stated to 
either the so called classical plateau state (CP, where 
the spin configuration resembles an Ising up-up-down 
state ttl) or the quantum plateau state (QP, where the 
spin configuration resembles a quantum singlet-up state 
••I). Experimental and numerical evidence for a quan- 
tum plateau at M = 1/3 was recently presented^ for 
Cu^(P20qOH)2, a newly synthesized compound that is 
very well described by spin 5=1/2 antiferromagnetic 
chains with period three modulated exchange couplings. 

Some insight about the magnetoelastic ground state 
can be obtained from the mentioned fixed modulation 
results at M = 1/3. When one considers a J\ — J 2 
chain with only nearest neighbors spin-phonon coupling, 
a lattice deformation that brings closer two neighbors 
to the same site (see Fig. [TJ upper panel) induces a 
spin exchange modulation in J\ forming open trimers. 
The ground state of the isolated trimer with S z = 1/2 
indicates^ the pinning of one of the classical plateau 
states, namely that with |J,| order on each trimer. 
Instead, if two of every three sites group together form- 
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FIG. 1: The upper panel describes a lattice deformation that 
enhances trimers; the lower panel corresponds to dimer en- 
hancement. 

ing dimers (see Fig. [IJ lower panel), the chain is driven 
to the quantum plateau state, that with spin singlets at 
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each dimcr. A different situation arises when one con- 
siders aiso next-nearest neighbors spin-phonon coupling, 
leading to modulations of both J\ and J2 exchanges. The 
particular modulation discussed in Ref. 0] could be ob- 
tained (see again Fig. [1] upper panel) when J\ and J2 
are modified so as to form closed trimcrs. In contrast 
with the previous example, we have shown that in this 
case trimer enhancement drives this system towards a 
quantum plateau state. A natural question is then which 
magnetic configuration corresponds to a given lattice de- 
formation in the general case. 

Motivated by the preceding discussion, we investigate 
in this paper the possibility, suggested by the present 
authors and collaborators in Ref. [H, of a spin-Peierls 
like displacive transition in an antiferromagnetic S = 1/2 
magnetoelastic J\ — J2 Heisenberg chain, when magne- 
tization is set to M=l/3 by an external magnetic field. 
We explore such a system with both nearest neighbors 
(NN) and next-nearest neighbors (NNN) spin-phonon 
couplings in the adiabatic approximation, allowing for 
modulations of J\ and J2 exchanges. This approach fol- 
lows the recent discussion in Ref. |ll|, where antiferro- 



magnetic zig-zag spin chain compounds such as CuGeO^ 
and L1V2O5 are argued to present NNN spin-phonon in- 
teractions at least of the same order as the NN ones; 
indeed, a numerical study of such magnetoelastic zig-zag 
chains at zero magnetization has lead to novel tetramer- 
ized spin-Peierls like phases. 

We show that the magnetoelastic ground state at zero 
temperature indeed favors several period three distortion 
patterns, stemming from a competition between clastic 
energy loss and magnetic energy change. These patterns 
spontaneously break translation symmetry, with differ- 
ent phases depending on the frustration ratio J2/J1 and 
the value of spin-phonon couplings. As we discuss be- 
low, there are essentially four different situations that 
arise when a lattice deformation of period three gener- 
ates a spin exchange modulation at M =1/3: the lattice 
shows two kinds of period three deformation patterns, 
namely (i) one tending to group three consecutive lattice 
sites into trimers and (ii) another one tending to group 
two of every three sites into dimers (see Fig. [T|). For 
each deformation pattern, depending on the microscopic 
parameters, the spin sector adopts either (a) a classical 
plateau configuration, pinned in the lattice with the J, 
spins sitting in the most convenient sites, or (b) a quan- 
tum plateau state, with the spin singlets located at some 
convenient links. A rich phase diagram is built, includ- 
ing all of the combinations of dimcr-like and trimcr-likc 
deformations with both classical and quantum plateau 
states. 

The paper is organized as follows. In Section II we 
present the model and its analytical treatment. The 
spin sector is described within the bosonization approach, 
while the phonon sector is described in the adiabatic ap- 
proximation by classical static deformations. In Section 
III we analyze this effective description by considering 
all of the relevant perturbation terms as semiclassical po- 



tentials, and draw a qualitative phase diagram with our 
results. Special emphasis is put on the characterization 
of the ground state phases that result from the combina- 
tion of frustration and magnetoelastic effects in different 
parameter ranges. In Section IV we present the results 
of Lanczos exact diagonalization of small systems, sup- 
porting the bosonization results. Finally, in Section V we 
present a summary and conclusions of the present work. 



II. DESCRIPTION OF THE MODEL AND 
BOSONIZATION APPROACH 

We consider the lattice Hamiltonian of a frustrated 
spin 5=1/2 Heisenberg chain, which can be written as 



Hm — 



ri+2 S n .S„ + 2 
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where S n are spin operators at site n and J n ,n+a > 
are antiferromagnetic NN (a = 1) and NNN (a = 2) 
spin exchange couplings. An uniform external magnetic 
field is also coupled to the spins in order to produce a 
global magnetization M = 1/3 (AI = 1 corresponding to 
saturation) . 

The interaction of spins in a homogeneous zig-zag 
chain (J n , n +i = Ji, Jn,n+2 = J2) with phonons is usually 
modeled by a linear expansion of the exchange couplings 
around the non distorted values J\ and J2 

J n ,n+1 ~ JlO- — A(u n+ \—U n )), 

Jn,n+2 ~ J 2 (l- B(u n+2 -Un)), (2) 

where u n is a scalar relevant coordinate for the displace- 
ment of ion n from its equilibrium position, and A, B are 
called the spin-phonon couplings at NN and NNN sites. 
The total Hamiltonian, including the elastic energy in 
the adiabatic approximation, is written as 

H T = -K^2(u n+1 - u n ) 2 + 

n 

+ S„.S, 1+ i + J2 S„.S„ + 2| — 

n 

- ^2 j^ 1 A(u n+1 - u„)S„.S n+ i + 
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(3) 



where K is the homogeneous spring stiffness. The 
first line corresponds to classical phonons elastic energy 
(Hep), the second one to the homogeneous magnetic 
Hamiltonian (Hm) and the rest to the spin-phonon in- 
teraction (H]), 



Ht — Hi 



CP 



H 



H, 



(4) 



Dimcnsionless parameters, convenient for numerical anal- 
ysis, are used below. They arc introduced using J\ as 
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the energy scale as follows: A -> A = (A/Kf^A, 
B -► B = (.h/R'Y^B, u„ -> u n = (K/J 1 ) 1 / 2 u n and 
Ji, J 2 — > a = J 2 /Ji- 

In order to observe semi-quantitatively the low energy 
properties of the model given by Eq.([3]), we employ the 
bosonization method which is generally powerful for the 
description of one-dimensional spin chains (see for in- 
stance Ref. fl3p. 

We start with the homogeneous magnetic Hamiltonian 
Hm- To obtain the corresponding low-energy theory 
one first applies the exact Jordan- Wigner transformation 
mapping spins onto spinless lattice fermions ij) n , then one 
introduces a continuum coordinate x = na with a the 
lattice spacing and writes a linear approximation for the 
low energy degrees of freedom around the Fermi level in 
terms of left and right-moving continuum fermions; the 
Fermi wave- vector Rf depends on the magnetization. For 
M = 1/3 one gets kp = 7r/3a, then 

yj n « j™^ R {na) + er m l^ L {na). (5) 

The continuum fermions are spinless and masslcss, allow- 
ing for Abelian bosonization; the complete Hamiltonian 
is finally mapped into a Gaussian term 

\ J dx[^-(d<f>) 2 +K L (d^ 2 ] (6) 

plus several vertex operators that are kept only when 
they are commensurate (non-oscillating in space) and 
constitute relevant perturbations to the Gaussian confor- 
mal field theory. Here <j> is a compactificd boson field de- 
fined on a circle, (f> = cfr+y^, and <p is its dual field defined 
by d x 4> = dt<p. The parameters v and Kl (Fermi velocity 
and Tomonaga-Luttingcr parameter respectively) depend 
on the microscopic parameters of the lattice Hamiltonian 
Hm; v is proportional to aJi, while Kl is dimensionless. 

A particular feature of the AI = 1/3 situation is 
that kp = 7r/3a makes commensurate a triple Umklapp 
process^, providing a perturbation term of the form 

- I <fecos(3V4^>) (7) 

in Hm- The coefficient g% is non-universal and, as well 
as v and Kl, depends on the renormalization group pro- 
cedure. 

It is known numerically^^ that the homogeneous mag- 
netic Hamiltonian Hm describes a gaplcss Tomonaga- 
Luttinger (TL) phase for < J 2 /Ji < a c = 0.5621. For 
<^c < J2/J1 ;$ 1-25 there exists a strong magnetization 
plateau at M = 1/3. Comparison of bosonization with 
these results shows that the the Tomonaga-Luttinger pa- 
rameter should be Kl > 2/9 for J2/J1 < a c , as this ren- 
ders the perturbation in Eq. ([7]) irrelevant. Then the co- 
efficient <?3 flows to zero under the renormalization group 
and the effective theory describes a gapless TL phase. 
On the other hand, for a c < J2/J1 < 1.25, it should be 
Kl < 2/9, making Eq. ([7]) a relevant perturbation. Thus 



this term opens a magnetic gap and explains the mag- 
netization plateau^ observed in this range. Moreover, 
the plateau ground state is known to be three-fold de- 
generate, with translation symmetry spontaneously bro- 
ken to an up-up-down configuration^; such configurations 
are described by the pinning of the bosonic field in one 
of the minima of Eq. ([7]) considered as a scmiclassical 
potential energy^, provided that c/3 > 0. We will in 
consequence qualitatively represent the plateau by the 
behavior of the non-universal coefficient 173 > as being 
smooth and non- vanishing only for a c < J2/J1 < 1.25, 
with a maximum at some intermediate value of J2/J1. 
We will not study here the regime J2/J1 > 1.25, where 
the M — 1/3 plateau is not present; this should be bet- 
ter done by starting with two spin chains with strong 
exchange J2, weakly coupled by a zig-zag interaction J\. 

Next, we consider the lattice deformations. From the 
knowledge of the M = 1/3 plateau magnetic ground state 
in the homogeneous J\ ~ J 2 chain with J2/J1 > ol c , one 
can argue that an adiabatic lattice deformation caused 
by the spin-phonon coupling in Eq. ([3]) will have period 
three. This is also supported by bosonization, as such 
a deformation is commensurate with kp = it /3a. More- 
over, as discussed in Ref. even for J2/J1 < a c period 
three deformations cause commcnsurability of relevant 
perturbations at M = 1/3 and provide a mechanism for a 
spin gap (magnetization plateau) in this regime. Numer- 
ical evidence of the dominance of period three lattice de- 
formations, obtained from self consistent computations, 
was also given in Q . A uniform deformation, leading to 
global size change, can also appear—; this would produce 
a uniform shift in J\ and J2, which is unessential to our 
present analysis. 

We will consider in this paper the most general period 
three deformation, without collective displacement, given 
by 

«n = sin(— n-x), (8) 

with the amplitude uq and a relative phase x as free pa- 
rameters. Our purpose is to search for the deformation 
that minimizes the magnctoclastic energy. Once a mini- 
mum of the total energy is found, the amplitude uq will 
indicate the deformation strength and the phase \ will 
relate the deformation pattern to the corresponding spin 
ground state characterized by the value of 4> at the po- 
tential minimum. 

From Eq. (0) the distortion of the NN bond length 
between sites n and n + 1, denoted by 5 n = u„+i — u n , 
is parameterized by 

2tv 1 

5 n = Mo cos (— (n + -) -x), (9) 
while the NNN distortion is given by 

8 n+ i + 8 n = u cos (— (71 + 1) - x)- (10) 
The elastic energy cost associated to deformations in Eq. 
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reads simply 



Hcp/Ji = -Nul 



(11) 



Finally, we consider the spin-phonon interaction 
Hamiltonian Hi induced by lattice deformations in Eq. 
([8]). Following the bosonization procedure one generates 
an extra renormalization of v and Kl and perturbation 
terms of the form 

J dx (A cob(V^</> + x ) + h cos(2V4^> - x )) , 

(12) 

thus introducing first and second harmonics of the bo- 
son field with coefficients proportional to the deformation 
amplitude uq. Notice that these operators are more rel- 
evant than the third harmonic in Eq. ([7]). and should be 
kept as well in the J2/J1 > ct c regime as for J 2 /J\ < a c , 
as far as Kl < 1/2. Even though the coefficients are 
non-universal and subject to renormalization. it is use- 
ful to report that a first order perturbativc computation 
yields h ~ A{l-C x J 2 /Ji) and / 2 ~ -A(l + C 2 q J 2 /Ji), 
where q = B /A and C\ , C 2 are positive constants with 
C2 < Ci. For a qualitative description, we will assume 
that /1 and f 2 depend on the microscopic parameters as 
suggested by these bare expressions. 

Putting all together, we can write the complete effec- 
tive theory as 



Ht = H C p + Hf ree + V e ff 



(13) 



where Hep is the classical elastic contribution given in 
Eq. dHJ), 

Hfree = \ J dx W) 2 + K L (04>) 2 ] (14) 
is the Gaussian part of the compactificd boson action and 



Veff = 



dx [uq fi cos(v 47T0 + x) 



2ir 2 a 2 

+uq f 2 cos(2V 4-7T0 — x) — 93 cos(3V Aire, 



(15) 



is the bosonic self-interaction potential defining a triple 
sine-Gordon theory^. Extensive analysis of competition 
between harmonics in multi-frequency sine-Gordon the- 
ories has been performe d 16 ' 17 ^ 18 , mainly focused on the 
double sine-Gordon model. The three-frequency case has 
also been recently discussed 1 ^. For our purpose it will be 
enough to perform a semiclassical treatment, as detailed 
in the next section. 



III. SEMICLASSICAL ANALYSIS OF THE 
EFFECTIVE THEORY 

The aim of the present work is to search for the pos- 
sibility of elastic deformations that lower the magne- 
toelastic energy with respect to the homogeneous non- 
deformed case. The simplest analysis of the effective the- 
ory obtained in the previous section, which has proved 



to be useful in related case a 7 ' 8 ' 14 , consists in treating the 
self-interaction terms in Eq. (fl~5|) as a classical potential 
to be evaluated in constant field configurations. 

Within this approximation the energy per site depends 
on three configuration parameters, uq, (f> and Xi an d is 
readily evaluated to 

e(uo, <t>, x) = -7^7 = 7«0 " |^ cos(3V4^,) (16) 



Ji N 



no 



+ 2^(/i cos(V4^> + X) + h cos(2V4^ - X )), 

so that the minima can be found analytically. Notice 
that this expression is invariant under simultaneous shifts 
V^7r0 — * V^ncj) + 27r/3, x ~^ X ~ 27r/3, in relation with 
the three equivalent locations of period three structures 
on the chain. This allows to restrict the analysis to < 
y/4n(f) < 2n/3 without loss of generality. Also a shift x 
X + 7T is equivalent to changing the sign of uo, allowing 
to consider < \ < ir. We report results within these 
restricted ranges. 

Among several local minima of the potential, the semi- 
classical energy is always found in one of the following 
situations: 



1. V^4> = 2tt/3, x = Jr/3, u = (/i + / 2 )/tt 2 , where 
the energy is evaluated to 



33 (/i + h? 



2tt 2 



4tt 4 



(17) 



2. V4^0 = tt/3, x = 2tt/3, mo 
the energy is evaluated to 



(/i - /2)/tt 2 , where 



.93 (A ~ hf 
2tt 2 



4tt 4 



(18) 



Before drawing a phase diagram, we discuss the phys- 
ical content of the possible phases. Following the usual 
bosonization rules to map </> to spin variables 2 ^, the value 
y/Iircj) = 27r/3 in the first solution indicates that the 
spin sector adopts a classical plateau configuration CP—, 
which corresponds to selecting one of the If J, degener- 
ate ground states of the homogeneous chain plateau. We 
will call its energy ecp- The relative phase x = t/3 plays 
together with the sign of uo in determining the elastic de- 
formation. For fi + f 2 >0 one finds a trimer-like elastic 
deformation grouping blocks of three spins (T, see Fig. 
[TJ upper panel); in the opposite case a dimcr-likc defor- 
mation is set, alternating two closer spins with a more 
separated one {D, see Fig.[TJ lower panel). In the second 
solution, the value V47T0 = 7r/3 is not one of the minima 
of the homogeneous chain potential, but it signals that 
the spin sector adopts a state that enhances quantum 
singlets in a •• | quantum plateau configuration 7 - QP. 
The corresponding energy will be called cqp. The rela- 
tive phase x = 27r/3 in this solution indicates that the 
lattice deformation is dimer-like (D) for /i — f 2 > and 
trimer-like (T) otherwise. 
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Depending on the coefficients /i, f 2 and 173, which in 
turn depend on the microscopic parameters, one of these 
solutions is selected as the global minimum and deter- 
mines the magnetoelastic ground state phase. 

In order to present a schematic phase diagram, we as- 
sume the qualitative phenomenological dependence of /1 , 
f 2 and (73 on the microscopic parameters detailed in the 
previous section. Following Ref. (TTf| we have chosen a 
ratio B = 1.5 A, which is used in the rest of the paper, as 
representative of materials where the NNN spin-phonon 
coupling plays an important role. The magnetoelastic 
phases found are shown in Fig. [51 We have checked that, 
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FIG. 2: Schematic magnetoelastic phase diagram with NN 
and NNN spin-phonon couplings related by B = 1.5 A. The 
classical up-up-down and quantum plateau phases are labelled 
by CP and QP respectively. The dimer and trimer elastic 
phases are labelled by D and T. Magnetoelastic patterns in 
each phase are shown by pictorial diagrams. 

within our approximations, all phase transitions result 
from level crossing of the above described local minima 
and can be then classified as first order. 

For the sake of illustrating the analysis leading to 
Fig. [21 we show in Fig. [3] the evolution of the energies 
e CP, e QP as functions of the NN spin-phonon couplings 
A, B = 1.5 A, fixing J2/J1 = 0.7. This situation lies well 
inside the homogeneous plateau regime, J 2 / J\ > ct c . The 
coefficients fi , f 2 are evaluated according to the first or- 
der bare result given in the previous section. The level 
crossing at A c ~ 0.65 shows the transition from CP mag- 
netic phase to QP. We remark that this transition is 
very different than that recently observed by the authors 
in [U, where the system passes from CP to a Z 2 broken 
symmetry phase and only then to a QP phase through 
an Ising like second order transition. 

Within each magnetic phase one can also identify the 
different clastic phases. Given A, we find critical val- 
ues of J2/J1 where uq changes sign: in the CP phase, 
the equation /i(a c i) = — /2(a:ci) defines a critical line 
•h / Ji = oid such that for J2/J1 < a c \ the system adopts 
a trimer like lattice distortion T while for J 2 j ' J\ > a c i 
the deformation is dimer like D. In contrast, in the QP 
phase we find a critical line a c2 where fi(a c2 ) = /2(ck C 2), 
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FIG. 3: Semiclassical energies for the classical and quantum 
plateaux minima in terms of the spin-phonon coupling A, for 
B = 1.5 A and J2/J1 = 0.7. The inset shows the ground state 
energy obtained by exact diagonalization of a system with 
N = 24 sites, after fitting the zero-energy level. 

being the lattice distortion of type D for J 2 j Ji < a c2 and 
of type T for J2/J1 > ct c2 . Using the bare expressions 
for /1, f 2 , these critical lines do not depend on A. 

A similar analysis can be made for J2/J1 < ol c . The 
most important difference is that in this region the 
Tomonaga-Luttingcr parameter is Kl, > 2/9 and the 
third harmonic is irrelevant. We represent this situation 
by setting c/3 = 0. As mentioned before, the magneti- 
zation plateau at M = 1/3 is induced by the coupling 
to the lattice^ through the first and second harmonics in 
Eq. (fl5|) . Unlike the previous case, there is no level cross- 
ing between cqp and ecp', the absolute energy minimum 
always corresponds to £qp, selecting the QP magnetic 
phase. We also find (fi—fz) > in the whole region, so 
that the elastic phase is of type D. 

The relative position of the elastic deformations and 
the magnetization profile at each phase is determined 
by the corresponding values of <p and Xi as can be found 
using Eq. (jSJ) and bosonization formulae. The four phases 
are described by diagrams in Fig. [2] 

Different ratios B/A can be analyzed similarly. We 
have observed that lowering this ratio produces an in- 
crease in the region characterized by the classical plateau 
and trimcr-likc deformations, with higher values of both 
a c i and A c . 

It is important to notice that the deformation am- 
plitude uq is proportional to /1, f 2l which are in turn 
proportional to the dimcnsionlcss spin-phonon couplings 
A. B. In Fig. [21 the elastic pattern evolves to the ho- 
mogeneous limit as A — > 0. By construction, the ef- 
fective theory in Eqs. (fl"4l [T5| describes in this limit a 
Tomonaga-Luttinger phase for J 2 / J\ < a c and a gapped 
sine-Gordon phase with triple degenerate ground state 
for J2/J1 > a c . 
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It is also interesting to mention that if our approach 
remains valid in the limit J2/J1 — > 0, describing a single 
NN spin chain, the system adopts a dimer like elastic de- 
formation which in turn induces a trimerized modulation 
with one larger and two smaller NN spin exchanges (c.f. 
Fig. [TJ lower panel). This model was recently studied 
by quantum Monte Carlo simulations of large systems in 
connection with Cu3(P20qOH)2t^, finding exactly the 
quantum plateau magnetic phase predicted by our anal- 
ysis. 

IV. NUMERICAL ANALYSIS 

In order to support the bosonization results in the pre- 
vious section, we performed a numerical analysis of the 
Hamiltonian in Eq. ((3]) by exact diagonalization of small 
clusters of size up to N = 24 sites with periodic boundary 
conditions. 

The strategy is the following: period three elastic de- 
formations without collective displacement are parame- 
terized by two independent bond distortions, say 81 = 
1*2 — ui and 82 = M3 — "2, while £3 = —Si — S2 and 
8n+3 = 8 n . For given values of J2/J1 and A, fixing B/A 
and M = 1/3, we computed by Lanczos diagonalization^ 
the exact ground state energy of the total Hamiltonian in 
Eq. (J3J) for a wide range of elastic deformations (S\ , 82, #3) 
and then selected the absolute minimum. 

We found that, in accordance with bosonization re- 
sults, the lowest energy configuration is always obtained 
(except for equivalent lattice translations) at one of lat- 
tice distortions patterns shown in Fig. [TJ 

1. (81,82,83) = (— |A, — ^A, A) that corresponds to 
the trimer-like phase (T), or 

2. (81,82, S3) = ( I A, i A, —A) that corresponds to the 
dimer-likc (D). 

In order to characterize the magnetic phases, we also 
computed the local magnetization profile < 5* > for the 
ground state. The order parameter 

Ms = jjYl cos ( Y {n ~ 2)) < 5 " > (19) 

n 

introduced in Ref. which is positive for the quantum 
plateau (QP) configuration and negative for the classical 
plateau (CP), is used to report the results. We found 
that both magnetic phases are realized at some region of 
either the T or the D elastic phases. 

A thorough scanning of the A > 0, J2/J1 > plane 
was made for TV = 24 sites, keeping B = 1.5 A. 
The magnetoelastic phases found are shown in Fig. [4] 
Representative scans at J2/J1 = 0.3, 0.6, 0.9 are shown 
in Fig. [5j showing 81, 82, 83, and Ms as functions of 
A. Notice that the deformation amplitude decreases for 
small A (a limit that corresponds to large stiffness K); 
due to finite size effects^ there is no deformation below 
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FIG. 4: Magnetoelastic phase diagram obtained by exact di- 
agonalization of N — 24 sites. Spin-phonon couplings ratio is 
set to B/A = 1.5. 

some finite value of A, which is seen to diminish with size 
by comparing N = 12, 18, 24 sites. The J2/J1 = 0.3 scan 
shows the phase D, QP for all A. In the J2/J1 = 0.6 case 
one can clearly observe the first order transition at some 
value of A (which in general depends on J2/ J\), with fi- 
nite jumps both in the deformations and the magnetic 
order parameter, from the T, CP to the D, QP phase. 
The same happens in the J2/J1 = 0.9 case, with a tran- 
sition from the D, CP to the T, QP phase. The region 
0.8 < J2/J1 < 0.9 shows that the critical line for transi- 
tion between T and D phases slightly depends on J2/J1', 
comparison with Fig.[5]mdicatcs that rcnormalization ef- 
fects on the bare coefficients j\ , /2 are not strong enough 
to impede our qualitative bosonization analysis. 

We have also analyzed different spin-phonon couplings, 
confirming the bosonization prediction that lowering the 
ratio B/A produces an increase in the region character- 
ized by the classical plateau and trimer-like deformations 
(c.f. Fig. [2]), with higher values of both a c i and A c . 

In summary, the numerical results confirm the semi- 
classical analysis given in previous section. 

V. CONCLUSIONS 

In the present work we have shown the existence 
of different magnetoelastic phases in S = 1/2 zig-zag 
antifcrromagnctic Ji — J 2 chains coupled to adiabatic 
phonons through nearest and next-nearest neighbor spin 
exchanges, at M = 1/3 magnetization. At zero tem- 
perature this situation corresponds to a magnetization 
plateau, either existing for the non-distorted homoge- 
neous chain with high enough frustration^iii or induced 
by spin-phonon coupling at lower frustration^. 

We performed a semiclassical analysis of the bosonized 
effective theory, supported by numerical exact diagonal- 
ization of small clusters up to 24 spins. We found that 
several spin-Peicrls like phases describe the ground state 
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FIG. 5: Lattice distortions 5i,2,3 and order parameter Ms as 
a function of A for B/A — 1.5 and J2/J1 = 0.3 (upper panel), 
0.6 (middle panel), 0.9 (lower panel). Each panel corresponds 
to a vertical scan in Fig. f4] Small A is not accesible due to 
finite size effects. 



of the system, depending on the microscopic parameters 
J2/J1 and spin-phonon couplings. In each of these phases 
a non trivial elastic deformation is favored, grouping to- 
gether blocks of two or three spins, while the magnetic 
sector adopts classical or quantum plateau states. 

A detailed analysis of a particular case, chosen as rep- 
resentative of materials with large ratio of next-nearest 
to nearest neighbors spin-phonon couplings^, shows the 
following magnetoelastic phases at zero temperature: 

(i) an up-up-down magnetic phase with a trimer-like lat- 
tice distortion when frustration is just enough to produce 
the M = 1/3 magnetization plateau in the homogeneous 
chain and spin-phonon couplings are low. 

(ii) an up-up-down magnetic phase with a dimcr-like lat- 
tice distortion for low spin-phonon couplings and higher 
frustration. 

(iii) a quantum plateau magnetic phase with dimcr-like 
lattice distortion for large spin-phonon couplings and low 
frustration. This phase is present even for such low frus- 
trations that would not produce a magnetization plateau 
in absence of spin-phonon interaction. 

(iv) a quantum plateau magnetic phase with trimer-like 
lattice distortion for large spin-phonon couplings and 
high frustration. 

Once the existence of non trivial magnetoelastic phases 
at zero temperature is proved, a natural question is to 
analyze the possibility of a spin-Peierls like transition in 
three dimensional materials with quasi-one-dimensional 
magnetic structure. Since a high temperature phase is 
expected to recover translation invariance, such a tran- 
sition could take place at some finite temperature^ , 
while an external magnetic field maintains the magne- 
tization M = 1/3. However, a finite temperature study 
should also take into account the eventual smoothing of 
the magnetization plateau. The critical temperature and 
the behavior of thermodynamic functions at the conjec- 
tured transition is suggested for future investigation. 
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